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Abstract. We review the status of lattice calculations of the deep-inelastic 
structure functions of the nucleon. In addition, we present some results on the 
pion and rho structure functions. 


INTRODUCTION 

The calculation of the deep-inelastic structure functions of hadrons from 
QCD is basically a non-perturbative problem. Perturbation theory can only 
describe the evolution of the structure functions from one scale to the other, 
provided one is at high enough scales. The most promising tool to solve this 
problem is lattice gauge theory. 

Three years ago we have initiated a program to compute the lower moments 
of the structure functions of the nucleon on the lattice. First results for the 
unpolarized structure functions F x , F 2 and the polarized structure functions g x , 
g 2 are now available [1-3]. It has been found that with the help of dedicated 
computers one is able to compute the lower moments of all four structure 
functions with satisfactory precision - at least for the non-singlet leading twist 
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contributions which we have addressed so far. It is needless to say that all these 
calculations have been performed in the quenched approximation, where the 
effect of internal quark loops has been neglected. Three or four moments are 
generally sufficient to determine the individual quark distribution functions [4]. 

Recently, we have extended our calculations to the structure functions of 
the pion and the rho [5]. The pion structure function is being measured at 
HERA, and data should soon become available. The structure functions of 
the rho are interesting because they give qualitatively new information about 
quark binding effects. 

Lattice calculations are subject to systematic errors arising from the finite 
lattice spacing, the finite volume, the extrapolation to the chiral limit and, of 
course, quenching. Before we can trust our results entirely, it is important to 
examine these effects carefully. This will be a major task in the future. 

The aim of this talk is twofold. In the first part we will present and discuss 
our results obtained over the last two years. This will be mainly about the 
nucleon structure functions and briefly about the structure functions of the 
pion and the rho. In the second part of the talk we will report on first attempts 
to reduce discretization errors. 


BASICS 


The theoretical basis of the calculation is the operator product expansion. 
For large Q 2 the moments of the nucleon structure functions are given by 
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where c\, c 2 and e\, e 2 are the Wilson coefficients, and v n , a n and d n are 
forward nucleon matrix elements of certain local operators O. For details see 
(e.g.) ref. [1], In parton model language 


Vnh = (x n ) U) = d xx n ^[q\ f \x) + q{ f \x)}, f = u,d, 
Vn+i = {x n ) {9) = [ d xx n g(x), 
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with A= A q(f\ and q-\(x), qi(x) being the quark distribution functions 
with spin f, j. with respect to the direction of motion. In eqs. (2), (4) we have 
omitted the sea quark contributions because we are working in the quenched 
approximation. The matrix elements d n have twist three and so have no parton 
model interpretation. 

For a spin-zero target like the pion one has q\(x) = qi(x). For a spin- 
one target like the rho the quark-spin distribution will depend on the spin 
projection m of the particle. Writing q m (x ) = 1/2 [qT(x) + q™(x)], one finds 
the new structure function [6] 

bi(x) = q°(x) — q 1 {x). (5) 

On the lattice one computes bare operators O(a), where a is the lattice 
spacing. These operators are in general divergent and must be renormalized. 
One may define finite operators renormalized at the finite scale // by 

0{ n) = Z 0 {an,g{a))0{a). ( 6 ) 

The renormalization condition here must match the renormalization condition 
used to calculate the Wilson coefficients, so that the /^-dependence drops out of 
the product of Wilson coefficients and Zq' s. In principle the renormalization 
constants can be calculated in perturbation theory. The problem is only that 
lattice perturbation theory converges very slowly [7]. 

NUCLEON STRUCTURE FUNCTIONS 

The lattice calculation of structure functions divides into two parts. Part 
one is the calculation of the hadron matrix elements. Part two is the deter¬ 
mination of the renormalization constants Zq. Altogether we have considered 
0(20) different operators so far. On the hypercubic lattice one is limited to 
operators with spin < 4. In the notation of eqs. (2)-(4) this means n < 3. 

For nearly all of the operators the renormalization constants have been com¬ 
puted perturbatively to one loop order [8-10]. A non-perturbative calculation 
has only been started recently [11,12]. In the few cases where we can compare 
perturbative and non-perturbative results we see differences of not more than 
5% if boosted perturbation theory [7] is applied. 

The numerical calculations of the matrix elements have so far been done at 
one value of the coupling, /3 = 6 /g 2 — 6.0. At this coupling the lattice spacing 
is a ~ 0.1 fm. We did calculations on two volumes, 16 3 x 32 and 24 3 x 32, and 
for several quark masses ranging from « 30 to 200 MeV. It is important to 



TABLE 1. Results for the unpolarized structure func¬ 
tions in the chiral limit. The lattice numbers are com¬ 
pared with the CTEQ3M [13] fit of the valence quark 
distribution functions. For (x), which has been calcu¬ 
lated using different procedures, we have averaged over 
the results in ref. [2]. 


Moment 

Lattice 

(quenched, g 2 ss 5 GeV 2 ) 

Experiment 
(g 2 = AGeV 2 ) 

(*)(“) 

0.410(34) 

0.284 

(x)W 

0.180(16) 

0.102 

(x)^-(x)^ 

0.230(38) 

0.182 

{x 2 )W 

0.108(16) 

0.083 

{x 2 )^ 

0.036(8) 

0.025 

{x 3 )W 

0.020(10) 

0.032 

(x 3 )W 

0.000(6) 

0.008 

(x)^ 

0.53(23) 

0.441 


TABLE 2. Results for the polarized structure functions in the chiral limit. 
The lattice numbers are compared with experiment [14,15] and the LO [16] 
fit to the valence quark distribution functions. The latter are referred to 
as valence. The data for A^u and Arefer to ^ i 2 = lOGeV 2 . 


Moment 

Lattice 

(quenched, g 2 « 5GeT 2 ) 

Experiment 
(g 2 = 3 - 5 GeV 2 ) 

A u 

0.841(52) 

0.823 valence 

Ad 

-0.245(15) 

-0.303 valence 

9A 

1.086(67) 

1.26 

fo dx (9i ~ 9i) 

0.176(14) 

0.163(10)(16) 

A^u 

0.198(8) 

0.169(18)(12) 

A Wd 

-0.048(3) 

-0.055(27)(11) 

fo dxx 2 g[ p) 

0.0150(32) 

0.012(1) 

fo 1 d xx 2 g[ n) 

-0.0012(20) 

-0.004(3) 

fo' da :x 2 g ( 2 p) 

-0.0261(38) 

-0.006(2) 

fo dxx 2 g < f > 

-0.0004(22) 

0.005(8) 
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FIGURE 1. The it-quark distribution function extracted from the lattice data compared 
with the experimental valence quark distribution function. The figure is taken from ref. [4]. 


have data at many quark masses to do a reliable extrapolation to the chiral 
limit. The calculation of nucleon matrix elements, in particular of higher 
derivative operators, requires high statistics. Our calculations on the 16 3 x 32 
lattice involved 0(1000) gauge field configurations, and 0(100) configurations 
on the 24 3 x 32 lattice. The calculations on the large volume were nearly as 
expensive in computer time as the high statistics calculations on the smaller 
volume as we went to smaller quark masses. For the calculation of the gluon 
distribution function we even used 0(5000) configurations. By comparing the 
results on the two volumes, we find that our calculations on the 16 3 x 32 lattice 
do not suffer from finite size effects. All results in this and the next section 
are for Wilson fermions. 

Our results for the unpolarized structure functions arc listed in table 1, and 
in table 2 we give the results for the polarized structure functions. Recent 
work on A q by two other groups 2 [17,18] give similar numbers to ours. We 
compare the lattice results directly with experiment, where possible, and oth¬ 
erwise with the phenomenological valence quark distribution functions, namely 
CTEQ3M [13] for the unpolarized structure functions and LO [16] for the po¬ 
larized structure functions. 

Let us now discuss the results in some detail. We begin with the unpolarized 

2 ) These groups also compute the sea quark contribution which one might object to in a 
quenched calculation. 





structure functions in table 1. We see that the lattice prediction for the 
lowest moment, (x), is significantly larger than the lowest moment of the 
phenomenological valence quark distribution function, both for the u and the 
d quark. When one goes to the higher moments the agreement improves. 

What does this mean for the structure functions themselves? In [4] the 
moments from table 1 have been converted to ‘real’ structure functions. In fig. 
1 the result is shown for the u-quark distribution function. This is compared 
with the phenomenological curve. We see that the ‘lattice’ structure function 
is quite different in shape from the experimental one. The difference shows 
mainly at x values around 1/2. We do not expect that the missing internal 
qnark loops will have an effect at such large values of x. If the effect is real, 
it could mean that the experimental structure functions contain substantial 
higher-twist contributions. One would expect higher-twist contributions to 
show up predominantly at large x. 

The lowest moment of the gluon distribution function, ( x)^ 9 \ turns out to 
be consistent with experiment. The relatively large error has to do with the 
fact that the calculation requires a delicate subtraction between two terms 
similar in magnitude [3]. Quark and gluon contributions to (x) must add up 
to one, as a result of energy-momentum conservation. We obtain 

(x} iu) + (x)( d) + (x) i9) = 1.12(23). (7) 

Let us now turn to the polarized structure functions in table 2. Here we 
have a lot of experimental data to compare with, so that we do not have to 
rely only on the phenomenological distribution functions. A recent, and yet 
unpublished, piece of information comes from the SMC measurement of the 
charge asymmetry of fast pions in polarized muon-nucleon collisions [15]. In 
the asymmetry the effect of sea quarks drops out, and one is lead directly to 
the moments of the valence quark distribution functions 

A (1 >«= ^ai u) , A«d= ^af 0 . (8) 

The lattice calculation of and Ais reported here for the first time [19]. 
For the higher moments sea quark effects should not play any significant role 
any more because they are restricted to the region of small x. 

With two exceptions, the moments of the polarized structure functions agree 
surprisingly well with experiment. One exception is the axial vector coupling 
constant qa = A u — Ad. Experimentally, this is one of the best known quan¬ 
tities, and the sea quark contribution is expected to cancel out. The lattice 
result turns out to be quite a bit lower than the experimental value. The other 
exception is the twist-three matrix element d 2 , contributing to the second mo¬ 
ment of the structure function g 2 . Comparing existing data on rq and g 2 , one 
is lead to the result d 2 ~ 0, both for the u and the d quark. We find that 
d 2 ~ — a 2 in the chiral limit. In the heavy quark limit, on the other hand, d 2 



vanishes as one would expect. We have no explanation for this discrepancy. 
We do not think that it can be explained by finite cut-off effects and the effect 
of quenching. 

PION AND RHO STRUCTURE FUNCTIONS 

The pion and rho structure functions have been computed from 0(500) 
gauge field configurations on the 16 3 x 32 lattice. As before, the coupling was 
taken to be f3 = 6.0. 

For the pion we obtain in the chiral limit 

(x) = 0.274(13), ( x 2 ) = 0.107(35), ( x 3 ) = 0.048(20). (9) 

If we compare this result with indirect information from the Drell-Yan cross 
section [20], we find the same picture as for the nucleon: (x) comes out to be 
significantly larger than the phenomenological value, while ( x 2 ) and ( x 3 ) are 
consistent with the data. 

The unpolarized rho structure function looks very similar to the pion struc¬ 
ture function. For the moment A q (the counterpart of Ain eq. (4) with 
m — 1 instead of 1/2) of the polarized structure function we find 0.59(5) in¬ 
dicating that the valence quarks carry 60% of the spin of the rho. This is 
about the same fraction as one finds for the nucleon. The lowest moment of 
the polarized structure function b\ turns out to be positive and unexpectedly 
large, albeit with large statistical errors. A possible interpretation of this is 
that the valence quarks have a significant orbital angular momentum. 

RECENT DEVELOPMENTS 

Wilson fermions (which we are using here) give rise to systematic errors of 
O(a), while the standard gluonic action induces errors of 0(a 2 ) only. Therefore 
it is the fermionic action which is most in need of attention. From hadron mass 
calculations we know that discretization errors can have a significant effect on 
the results at present values of the coupling [21]. 

Since it is very expensive to reduce cut-off effects by reducing a, a better 
way of reducing them is by improving the action. A systematic improvement 
program reducing the cut-off errors order by order in a has been proposed by 
Symanzik [22] and was further developed in [23]. By adding the (irrelevant) 
term 


( 10 ) 

X 

to the (Wilson) action [24] and choosing the coefficient c appropriately [25], 
one can reduce the cut-off errors from 0(a) to 0(a 2 ), provided the operators 
are improved as well. 
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FIGURE 2. The moment (x) for u and d quark distributions as a function of the quark 
mass at (3 = 6.0. We compare results using the Wilson action (o) with results from the 
improved action (x). The lines are a linear extrapolation of the Wilson action data to the 
chiral limit. The extrapolated values (•) are compared with the CTEQ results (^). 
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FIGURE 3. The same as fig. 1 but for A q. Here the extrapolated values are compared 
with the phenomenological distribution LO of ref. [16]. 






In figs. 2 and 3 we show first results of the 0(a) improved theory for the 
moments 3 (x) and A q. The calculations are done on the 16 3 x 32 lattice at 
/3 = 6.0. We compare the improved results with the unimproved numbers. To 
facilitate the comparison, we have only looked at one operator on one lattice. 
In the case of A q we see no difference between the two results at any value 
of the (bare) quark mass. In the case of the moment (x) we see only a small 
effect, which is however not statistically significant. For both moments the 
quark mass dependence is rather weak and unspectacular, so that a linear 
extrapolation to the chiral limit is justified. 

We may conclude that our numbers for (x) and A q are not far from the 
continuum values (of the quenched theory). This means that we have a real 
problem with the unpolarized structure functions. It is very likely that the 
solution lies in unexpectedly large higher-twist contributions at our values 
of Q 2 . The situation of the axial vector coupling constant qa has also not 
changed. One should, however, keep in mind that the naive parton model 
identification of q a with A u — Ad might fail. 

Whether discretization errors are also small for the other moments remains 
to be seen. 


CONCLUSIONS 

We have reviewed the current status of lattice structure function calcula¬ 
tions. The results so far are encouraging. But it is clear that much more 
remains to be done. We have only begun to improve the calculation by sys¬ 
tematically removing discretization errors, and this will certainly occupy us 
throughout the next year. 

A topic which is of upcoming interest [29] is power corrections (see also 
above). This includes higher twist effects on the one hand and renormalon 
contributions on the other. We hope to be able to present first results on this 
topic at next year’s workshop. 
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